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Abstract. We study the symmetry properties of autonomous integrating factors from an 
algebraic point of view. The symmetries arc delineated for the resulting integrals treated as 
equations and symmetries of the integrals treated as functions or configurational invariants. 
£h ' The succession of terms (pattern) is noted. The general pattern for the solution symmetries 

for equations in the simplest form of maximal order is given and the properties of the 
associated integrals resulting from this analysis are given. 
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! 1 Introduction 

p . 

It is well-known that, when a symmetry is used to determine a first integral for a differential 
equation, the symmetry provides an integrating factor for the equation and remains as a sym- 
metry of the first integral. For first-order ordinary differential equations the direct determination 
, ^ ' of the integrating factor is known and algorithms for finding integrating factors for equations 

^ ■ of higher order have been developed. In 1999 Cheb-Terrab and Roche |5j presented a systema- 

tic algorithm for the construction of integrating factors for second-order ordinary differential 
equations and claimed that their algorithm gave integrating factors for equations which did not 
possess Lie point symmetries. In 2002 Leach and Bouquet [3] showed that for all equations 
except one of which Cheb-Terrab and Roche [2] had found integrating factors had symme- 
tries which were not necessarily point symmetries but generalised or nonlocal. In the same 
year, 2002, Abraham- Shrauner |3] also wrote a paper to demonstrate the reduction of order 
of nonlinear ordinary differential equations by a combination of first integrals and Lie group 
symmetries. The latter and former motivated us hereby to investigate the underlying proper- 
ties of autonomous integrating factors and the associated integrals treated as equations and as 
functions. Observations are made and inferred in general for any nth-order ordinary differential 
equation of maximal symmetry. These will also be extended to include other types of equations 
in a separate contribution. 

Program LIE |2j is used to compute the symmetries for the different cases considered. The 
knowledge of the symmetries of first integrals of the equation does give rise to some interesting 
properties of the equation itself. For example, the Ermakov-Pinney equation [7j which in its 
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simplest form is 

where K is a constant. In theoretical discussions the sign of the constant K is immaterial and 
in fact it is often rescaled to unity. 
The general form of (fT]). videlicet 

9 / \ 1 

occurs in the study of the time-dependent linear oscillator, be it the classical or the quantal 
problem, as the differential equation which determines the time-dependent rescaling of the space 
variable and the definition of 'new time'. Some of the references for this are [H1IH]. 

Another origin of Q — of particular interest in this work — is as an integral of the third-order 
equation of maximal symmetry which in its elemental form is y'" = 0. 



2 Equations of maximal symmetry 



Definition 1. We define a first integral / for an equation of maximal symmetry, E = y^ = 0, 
as I = f(y,y',y',.. . ,y (n_1) ), where 



dl 

dec 



E=0 



dx 



E=0 



This means that, if g{x, y, y' , y" , . . . , y^ n ^) is an integrating factor, then 



dl 

dx 



E=0 



gE (x,y,y',...,y 



(n) 



E=0 



0. 



We start by considering the well-known third-order ordinary differential equation of maximal 
symmetry 



III n 

y = 



(2) 



which has seven Lie point symmetries. These are 



Gi = d y , G 2 = xd y , G 3 = x 2 d y , G 4 = yd y , 
G$ = d x , G e = xd x + yd y , G 7 = x 2 d x + 2xyd y . 



(3) 



The algebra is {A\ © s sl(2, R)} © s 3A±. The autonomous integrating factors for (j2J are y" and y. 
We list the symmetries and algebra when each of the integrals is treated as an equation and as 
a function. 

When we multiply y'" = by the integrating factor y" , we obtain y"y"' = 0. Integration of 
this expression gives \ (y") 2 = k, where k is a constant of integration. This gives rise to three 
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cases which we list as follows: 

y" = o 

G\ = d y , 
G 2 = xd y , 

G 3 = yd y , 

Ga = d x , 
G 5 = xd x , 

G 6 = x 2 d x + xydy, 

G 7 = yd x , 

G 8 = xyd x + y 2 d y , 

and, when y" = k is treated as a function, we have 

Gi=d y , G 2 = xd y , G 3 = d x , G 4 = xd x + 2yd y . 

Remark 1. When y" = k is treated as an equation, we have two cases, that is, y" = and 
y" = k for which the algebra is sZ(3,R) : 2Ai S {s/(2,R) ®A 1 }®2A 1 Q3JQUIT3]. If y" = k is 
treated function, the algebra is A\ 9 : A 2 S 2 Ax jTTl fT2l fT8l fll] . 

If y is used an the integrating factor, we obtain yy'" = 0. Integration of this equation gives 
yy" ~ \y' 2 = k which can be written as (y 1 ^ 2 )" = k/(y 1 / 2 ^ and is the simplest form of the 
Ermakov-Pinney equation jSld- As before we write down the point symmetries corresponding 
to the three cases of the differential equation u" = k/u 3 , where u = y 1 ^ 2 . Program LIE |3] gives 



the following: 










u" = 




u" = k/u 3 




u" = k/u 3 


G\ = d u , 


Gi 


= d x , 


Gi 


= d x , 


G 2 = xd u , 


G 2 


= 2xd x + ud u , 


G 2 


= 2xd x + ud u , 


G 3 = ud u , 


G 3 


= x 2 d x + xud u , 


G 3 


= x 2 d x + xud u , 


Ga = d x , 










G 5 = xd x , 










Gq = x 2 d x + xud u , 










G 7 = ud x , 










Gs = xud x + u 2 d u - 











The transformation of yy" — ^y' 2 = k to u" = k/u 3 does not make a difference in terms of the 
symmetries as we just have a point transformation in this case. The other obvious integrating 
factors for @ are 1, x and \x 2 which give 

y", 

xy" - y', 
-x 2 y" — xy' + y. 

(Note that the numbering of the fundamental first integrals follows the convention given in 
Flessas et aZ H3H5I0 



y" = k 



G l = d, l 

G 2 
G 3 



Ga 
G 5 

G 6 
G 7 
Gs 



(^x 2 k - y^j d y , 



x 2 d x + ( xy + -x 3 k ) d y , 



y - -x 2 k ) d x - x 6 k 2 d y , 



1 



xy 



x 3 k ) d x + [y 



c A k 2 d 



iY" = o — h = 

x-y"' = — »• h = 
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• The integration of equation © , which is a feature of the calculation of the symmetries of all 
linear ordinary differential equations of maximal symmetry ^Uj , by means of an integrating 
factor gives a variety of results depending upon the integrating factor used. 

• The characteristic feature of the Ermakov-Pinney equation is that it possesses the three- 
element algebra of Lie point symmetries, sZ(2,R), which in itself is characteristic of all 
scalar ordinary differential equations of maximal symmetry. 

The fourth-order ordinary differential equation y lv = has autonomous integrating factors y' 
and y'" . If we use y' as an integrating factor in the original equation and integrate, we obtain 

y'y'" - \{v"? = *■ (4) 

Equation @ is a generalised Kummer-Schwartz equation for k = and for / a variation 
on the Ermakov-Pinney equation as it can be written in the form 



(yr /2 )"=(k/((y'f /2 )). 



The three cases for the integral in Q treated as an equation and as a function give the following 
results: 

y'y'" - \{y"? = o y'y'" - \{y"? = k y'y'" - \(y") 2 = k 

G\ = d x , G\ = d x , G\ = d x , 

G*2 = xd x , G2 = d y , G2 = d y , 

Gs = ydy, G 3 = xd x + 2yd y , G 3 = xd x + 2yd y , 

G4 = d y . 

The use of y'" as an integrating factor gives y'" = k. If k = 0, then we just have seven point 
symmetries as those of equation (f2"J). The two remaining cases give 

y'" = k y'" = k 

G\ = d x , G\ = d x , 

G 2 = xdy, G 2 = d y , 

G3 = ]jX 2 d y , G 3 = xdy, 

G4 = dy, G4 = X^dy, 



G 5 = xd x + -x 3 kd y G 5 = xd x + 3yd y , 



Ge= (y- ^ x3 kj d y , 
G 7 = x 2 d x + ( 2xy + \x^k \ d, 
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Consider the fifth-order equation of maximal symmetry given by 

y v = (5) 

with autonomous integrating factors y, y" and y" . If we multiply © by the first integrating 
factor and integrate, we obtain the integral 

yy"-y'y'" + \{y"? = k. (6) 
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We consider the three cases for @ treated as an equation with k = 0, k ^ and as a function. 
y y iv _ y , y m + 1 {y/l) 2 = yy iv _ yly ,n + 1 (y // )2 = fe yy iv _ + 1 (y „ )2 = fe 

Gi = d x , G\ = d x , G\ = d x , 

G 2 = xd x , G 2 = xd x + 2yd y , G 2 = xd x + 2yd y , 

G3 = yd y , G 3 = x 2 d x + Axydy, G 3 = x 2 d x + 4xyd y , 

G4 = x 2 d x + 4xyd y . 

Remark 2. For easier closure of the algebra in the first case xdx can be written as xd x + 2yd y . 



• We also observe that there is no difference when the integral is treated as a function and 
as an equation. It is important to note that, if y is an integrating factor of y^ = 0, then 
the integral obtained using this integrating factor always has the sl(2,~R) subalgebra. 

• We further observe that for the peculiar value of the constant, that is, k = 0, there is the 
splitting of the self-similarity symmetry into two homogeneity symmetries. 

The integrating factor y" with © gives the following results 

y"y iv - \{y"'? = o y"y iv - \{y"') 2 = k y"y iv - \{y'") 2 = k 

G\ = d x , G\ = d x , G\ = d x , 

G 2 = xd x , G 2 = d y , G 2 = d y , 

G3 = 8y, G3 = X8y, G3 = X8y, 

G4 = xd y , G4 = xd x + 3yd y , G4 = xd x + 3yd y , 

G 5 = ydy. 

If we use y lv as the integrating factor of © and integrate, we obtain 

y [v = k. 

We delineate the three cases below: 

y iv = y iv = k y iv = k 

G\ = d y , G\ = dy G\ = dy, 

G 2 = xdy, G 2 = xdy G 2 = xdy, 

G 3 = x 2 d y , G 3 = ^x 2 d y G 3 = x 2 d y , 

G 5 = ydy, G 5 =d x G 5 = d x , 

G 6 = d x , Gq = 6xd x + x 4 kd y G 6 = xd x + 4yd y , 

G 7 = xd x , G 7 = (24y-x 3 k)d y , 

G 8 = x 2 d x + 3xyd y , G 8 = 2Ax 2 d x + (72xy + x 5 k)d y . 

The differential equation 

y vi = (7) 

has integrating factors y' , y'" and y v . If we use y' as the integrating factor, we obtain 

/v // iv 1 / ///\2 r 

yy -y y +^{y ) =k 
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which leads to the cases below. 

y'y v -y"y [v + l(y"') 2 = o 

G\ = d y , 

G 2 = yd y , 
G3 = d x , 
G4 = xd x . 



y'y v 

Gi = 
G 2 = 
G 3 = 







d x , 

xd x + 3yd y , 



y'y v 

Gi = 
G 2 = 
G 3 = 



y"y w + - 2 {y'"f 



d 



d x , 

xd x + 3yd y , 



The use of y'" as the integrating factor for (J7J) leads to 
k. 



Ill V 

y y 



W v ) 2 



The three cases give the following results: 



y"Y 




y"'y v 


-\(y iv ) 2 = k 


y"'y v 


~\^) 2 = 


Gi = 


3y, 


G x = 


dy, 


Gi = 


dy, 


G 2 = 


d x , 


G 2 = 


d x , 


G 2 = 


d x , 


G 3 = 


xQy % 


G 3 = 


Xdy^ 


G 3 = 


Xdy-, 


G4 = 


X 3y , 


G 4 = 


X dy , 


G 4 = 


X dy , 


G 5 = 


ydy, 


G 5 = 


xd x + kydy, 


G 5 = 


xd x + kydy, 


Ge = 


xd x - 











If we use y v as an integrating factor, we obtain 

y v = k. 

We also have the three cases as mentioned above to obtain 

y v = k 

dy, 
d x , 

Xdy , 

x 2 d, 
x 3 d. 





y v = o 




G x 


= dy, 


G L 


G 2 


= ydy, 


G 2 


G 3 


= d x , 


G 3 


G 4 


Xdy^ 


G 4 


G 5 


X Oyi 


G 5 


Gq 


= X 3 dy, 


G 6 


G 7 


X dy. 


G 7 


G 8 


xd x , 


G 8 


Gg 


= x 2 d x + 4:xyd y , 


Gg 



X 



!)■■ 



xd x ~\~ Q^kx dy, 



Gi 
G 2 
G 3 
G 4 
G 5 
G 6 

G 7 



y ' 
d x , 

Xdy , 

x 2 d, 

3 a 



X 



.'/- 



X dy , 

xd x + 5yd y , 



kx 5 I d v , 

120 J y 

x 2 d x + (^xy + r^kx 6 ^ 



kx 6 I dy. 



For the differential equation y vn = we have the integrating factors y, y" , y lv and y vl . The 
integrals corresponding to these integrating factors respectively are 



yy 



y"y vi 



y lv y vl 



y'y v + y"y lv 



\{y'"? 



k, 



y'"y v + \(y iv ) 2 



(8) 



k, 



k. 
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If y is used as the integrating factor, we have the integral yy vl — y'y v + y"y" — \{y"') 2 
is treated as an equation for k = 0, k ^ and as a function. This gives the following results: 

G\ = d x , G\ = d x , G\ = d y , 

G 2 = xd x , G 2 = xd x + Sydy, G 2 = xd x + 3yd y , 

G 3 = yd y , G 3 = x 2 d x + 6xyd y , G 3 = x 2 d x + 6xyd y , 

G 4 = x 2 d x + Qxydy. 

The integral corresponding to the integrating factor y" leads to the following cases: 



k which 



y"y vi 


-y"'y v - 




1 = y"y v[ 


-y"Y + \{y [v Y 


1 = k y"y vi 


-y'"y" + \(y™) 2 


Gi = 


8y, 




G l = 


dy, 


Gi = 


By, 


G 2 = 


dx, 




G 2 = 


9 X , 


G 2 = 


d x , 


G 3 = 


Xdy, 




G 3 = 


Xdy , 


G 3 = 


Xdy , 


G 4 = 


xd x 1 




G 4 = 


xd x + 4yd y , 


G 4 = 


xd x + 4yd y , 



k 



G 5 = ydy. 

For the integrating factor y lv we have the cases: 



yiV yVi 






-\(y v ) 2 = k 


y iv y vi 




Gi = 


dy, 


G a = 


dy, 


Gi = 




G 2 = 


d x , 


G 2 = 


d x , 


G 2 = 


d x , 


G 3 = 


ydy, 


G 3 = 


Xdy , 


G 3 = 


Xdy, 


G 4 = 


Xdy, 


G 4 = 


X dy. 


G 4 = 


X dy, 


G 5 = 


X dy. 


G 5 = 


X 3 dy, 


G 5 = 


X 3 dy, 


G% = 


X 3 dy, 


Gq = 


xd x + 5yd y , 


G 6 = 


xd x + 5yd y , 


G 7 = 


xd x . 











The last of the four integrating factors y vl leads to y^ 
following results: 





y vl = o 




Gi = 


dy, 


<n = 


G 2 = 


d x , 


G 2 - 


G 3 = 


ydy, 


G 3 : 


G 4 = 


Xdy, 


G 4 : 


G 5 = 


X dy , 


G 5 - 


G 6 = 


X dy, 


G 6 : 


G 7 = 


X dy , 


G 7 : 


Gs = 


X dy, 


G 8 : 


G9 = 


xd x , 


G 9 : 


G10 = 


= x 2 + bxyd y , 


G10 



y 



k 



k. We have for the three cases the 
y vi = k 



d, 
d x 



xd, 

~^ 2 dy, 
^dy, 

2A x49y ' 

x d-r ~\~ , k x d>i 1 , 

120 y 



Gi 
G 2 
G 3 

G 4 
G 5 
G 6 
G 7 
Gs 



dy, 

d x , 

Xdy , 

x 4 <9. 



x 



5 d 



!)■■ 



xd x + 6yd y , 



V 
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The differential equation 

y viii = o (9) 

has integrating factors y' , y'", y v and y vn . If we use y' in @ and integrate the resulting equation, 
we obtain the integral 

y'y yii -y"y vi + y"'y v -\{y 1Y f = k. (10) 

The three cases of the integral in (|l(Jj) being treated as an equation with k = and k ^ 
and as a function are given respectively below: 



Gi = 


8y 


Gt=dy 


G\ = d y , 


G 2 = 


yd y , 


G2 = d x , 


G2 = d x , 


G 3 = 


d x , 


G 3 = xd x + 4yd y , 


G 3 = xd x 


G 4 = 


xd x - 






y'" is 


used as an integrating factor, 


we obtain 



.... 1 9 

/// Vll IV VI i / V \ z 7 

y y -y y +^{y ) =k 

with the following respective cases: 



Gt 


= 9y, 


Gi 


= 9y, 


Gi 


= 9y, 


G 2 


= ydy, 


G 2 


= d x , 


G 2 


= d x , 


G 3 


= d x , 


G 3 


Xdy) 


G 3 


Xdy , 


G 4 


— xdy, 


G 4 


X dy^ 


G 4 


— X dy-, 


G 5 


= X 2 dy. 











The use of y v as an integrating factor gives 

v°v*-\kry=k. (11) 

Equation (|11() is of the Ermakov-Pinney type. The three cases can be delineated as follows: 

k 



y v y vii 








y v y vii 




Gi = 


dy, 


Gi = 


dy, 


G 1 = 


dy, 


G 2 = 


ydy, 


G 2 = 


d x , 


G 2 = 


d x , 


G 3 = 


d x , 


G 3 = 


Xdy , 


G 3 = 


Xdy , 


G 4 = 


Xdy^ 


G 4 = 


X dy^ 


G 4 = 


X dy^ 


G 5 = 


X dy. 


G 5 = 


X dy^ 


G 5 = 


X dy. 


Gq = 


X Dy, 


Gq = 


X dy^ 


G 6 = 


X dy. 


G 7 = 


X dy. 


G 7 = 


xd x + 6yd y , 


G 7 = 


xd x + 6yd y 


G 8 = 


xd x . 
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If y v11 is used as an integrating factor in Q , we obtain y v11 = k with the following symmetries 
for each of the three cases: 

y vii = y vii = k y vii = k 

G\ = d y , G\ = d y , G\ = d y , 

G 2 = yd y , G 2 = d x , G 2 = d x , 

G 3 = d x , G 3 = xd y , G 3 = xd y , 

G^. = Xdyj G4. = ~^X 2 3y , G4. = X^dy, 

G5 = X^dy, G5 = — X^dy, G5 = X^dy, 

G 6 = x 3 d y , G 6 = 24 x/L dy, G & = x ^ y , 

G 7 = x^dy, G 7 = -^ x ' J dyi G 7 = x ' J d y , 

G 8 = x 5 d y , G 8 = ^ x %, G 8 = x%, 

G 9 = x 6 d y , G 9 = xd x + ^ kx7 dyi G 9 = x9 x + 7 vdy 



Gin = xd r , Gin = ( y kx 7 ) d,, 

iu x, ru yy 504Q j , 



Gn = x 2 d x + ( 6xy + -^^ k ) 

3 Relationship between fundamental integrals 
and integrals obtained from integrating factors 

Consider the example of the third-order ordinary differential equation y'" = with the three 
fundamental integrals together with the appropriate associated point symmetries from the sub- 
algebra sZ(2,R): 

G 7 = x 2 d x + 2xyd y , h = -x 2 y" - xy' + y, 

G 6 = xd x + yd y , h = xy" - y\ 

G 5 = d x , h = y". 

The numbering of the symmetries follows that of the listing of Lie point symmetries in © and 
the ordering of the integrals is in terms of their solution symmetries. Then the autonomous 
integral associated with the integrating factor y comes from the combination 

J = h I,- l -ll = yy"- l -y'\ 

Proposition 1. All the integrals obtained using y as an integrating factor always have the 
sl(2,M.) subalgebra whereas the fundamental integrals only have one of the sZ(2,R) elements. 

Proof. To prove the first proposition we consider the s/(2,K) subalgebra Ai = d x , A 2 = 
xd x + yd y and A3 = x 2 d x + 2xyd y and the fundamental integrals I\, I 2 and I 3 respectively. Then 
we have the following: 

Ai/i=/ 2) A 2 h = h, A 3 /i = 0, 
Ai/ 2 = I 3 , A 2 / 2 = 0, A 3 / 2 = -2/i, 

A1/3 = 0, A 2 / 3 = -7 3) A 3 7 3 = -2/ 2 . 
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We also observe that Aj J = for i = 1,2, 3. In fact it is easy to show that Aj J = eijkljh- 
This is shown below as follows: 

AiJ = hh - hh = 0, 
A 2 J = h h ~ hh = 0, 

A 3 J = -2I 2 h + 2hh = 0. ■ 
In general we have 

Im= £ ) ' . 1 „ y (n - J - 1) , i = 0,1,..., 71-1, (12) 

n — ? — i — 1 

so that for n = 3, /30 = ii, /31 = h and ^32 = I3. The symmetries Ai, A2 and A3 operating on 
the fundamental integrals then yield in general 

Ai-Zfii In,i+1: 
Mini = (1 - ~i)Ini, 

Mhi = -(n + i-3)(n- i)I n ,i-i, 

Inn — 0. 

If we take for example A^I^i = —i(3 — i)hi—i with n = 3 and i = 0, 1, 2, we obtain 



A3/30 = 


> 


A3J31 = - 


2I3O) 


A3 132 = - 


2%, 


where as above 


ho = h, hi = h and 7 32 = ^3 



Proposition 2 ([lf|). If we take the equation of maximal symmetry = 0, the sZ(2,R) 
subalgebra maps back to itself and is preserved. 

Proposition 3. For the fifth-order equation y v = the autonomous integral emanating from 
the integrating factor y can be obtained from J = IqI^ — I\I 3 + \l 2 , where 

j 1 4 iv ^ 3 /// 1 1 2 // / 1 

= 24 ^ V ~ 6 X V + 2 Xy ~ Xy y ' 

I 
h = \x 2 y" - xy'" + y", 
I 3 = xy iv -y">, 
h = y iv - 



Proposition 4. The fourth-order equation also has an autonomous integral J defined as J = 
hh — 2^2 > where 

h = \x 2 y"' - xy" + y', 
l2 = xy'"-y", 

h = y'"- 

Proposition 5. It can be shown that the differential equation y vl = also has the autonomous 
integral J which is defined as J = hh ~ hh + hh ~ \l 2 with the Ii (i = 0, 6) being redefined 
appropriately. 
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4 Conclusion 

If = f(x, y,y' , ■ ■ ■ , y n_1 ) is an nth-order ordinary differential equation and g(x, y,y' , . . . , 
yU-i)^ _ ^ j g an integral, the integral obtained by multiplying the equation by the integrating 
factor and integrating once possesses certain symmetries when treated as a function, an equation 
for the general constant and a configurational invariant (k=0). It is important to note that, if 
y is an integrating factor of y^ = 0, then the integral obtained using this integrating factor 
always has the s/(2,R) subalgebra whereas the fundamental integrals only have one of the 
sl(2,M.) elements. We further observe that for the peculiar value of the constant, k = 0, there 
is the splitting of the self-similarity symmetry into two homogeneity symmetries. The third- 
order ordinary differential equation is actually special and leads to the Ermakov-Pinney type 
equation. The fourth-order ordinary differential equation y lv = has y' as one of its autonomous 
integrating factors which leads together with the the original equation upon integration to 
the generalised Kummer-Schwartz equation. An extension to other types of equations will be 
completed in a separate contribution. The question of what Lie point symmetries of an ordinary 
differential equation are also shared by all its first integrals will form the basis for the next 
contribution. 
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